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Abstract 

The term active matter describes diverse systems, spanning macroscopic (e.g. shoals of fish 
and flocks of birds) to microscopic scales (e.g. migrating cells, motile bacteria and gels formed 
through the interaction of nanoscale molecular motors with cytoskeletal filaments within cells). 
Such systems are often idealizable in terms of collections of individual units, referred to as 
active particles or self-propelled particles, which take energy from an internal replenishable 
energy depot or ambient medium and transduce it into useful work performed on the environ- 
ment, in addition to dissipating a fraction of this energy into heat. These individual units may 
interact both directly as well as through disturbances propagated via the medium in which they 
are immersed. Active particles can exhibit remarkable collective behaviour as a consequence 
of these interactions, including non-equilibrium phase transitions between novel dynamical 
phases, large fluctuations violating expectations from the central limit theorem and substantial 
robustness against the disordering effects of thermal fluctuations. In this chapter, following a 
brief summary of experimental systems which may be classified as examples of active matter, 
I describe some of the principles which underlie the modeling of such systems. 

1 Introduction to Active Fluids 

Anyone who has admired the intricate dynamics of a group of birds in flight or the coordinated, 
almost balletic maneuvers of a school of swimming fish can appreciate the motivation for the 
study of "active matter": How do individual self-driven units, such as wildebeest, starlings, fish or 
bacteria, flock together, generating large-scale, spatiotemporally complex dynamical patterns [51 J? 
What are the rules which govern this dynamics and how do the principles of physics constrain the 
behaviour of each such unit? Finally, what are the simplest possible models for such behaviour 
and is there any commonality to the description of these varied problems Il46ll47ll48l ? 

Describing these diverse problems in terms of individual "agents" which evolve via a basic 
set of update rules while interacting with other agents provides a general way of approaching a 
large number of unrelated problems. These include the description of the propagation of infectious 
diseases in a population, the seasonal migration of animal populations, the collective motion and 
coordinated activities of groups of ants and bees and the swimming of shoals of fish [|48l . Of these 
problems, the subset of problems involving agents whose mechanical behaviour at a scale larger 
than the individual agent must be constrained by local conservation laws, such as the conservation 
of momentum, forms a special class p4|. It is these systems that are the primary focus of this 
chapter. 

Generalizing from the examples above, a tentative definition of active matter might be the 
following: Active matter is a term which describes a material ( either in the continuum or naturally 



decomposable into discrete units), which is driven out of equilibrium through the transduction 
of energy derived from an internal energy depot or ambient medium into work performed on the 
environment. Such systems are generically capable of emergent behaviour at large scales. The 
hydrodynamic description of such problems should be equally applicable to situations where the 
granularity of the constituent units is not resolved, but which are rendered non-equilibrium in a 
qualitatively similar way ||48ll39l . 

What differentiates active systems from other classes of driven systems (sheared fluids, sedi- 
menting colloids, driven vortex lattices etc.) is that the energy input is internal to the medium (i.e. 
located on each unit) and does not act at the boundaries or via external fields. Further, the direction 
in which each unit moves is dictated by the state of the particle and not by the direction imposed 
by an external field. 




Figure 1: A school of fish, illustrating the tendency towards parallel alignment while swimming. 
Picture courtesy Prof. R. Kent Wenger 

Such individual units are, in general, anisotropic, as in Fig.Q] Collections of such units are thus 
capable of exhibiting orientationally ordered states. A canonical example of an ordered state of 
orientable units obtained in thermal equilibrium is the nematic liquid crystal, in which anisotropic 
particles align along a common axis. Active systems of such units can, in addition, exhibit directed 
motion along this axis; a concrete example is illustrated in Fig. [TJ In the context of active particles 
in a fluid, the terminology "swimmers" or "self-propelled particles" is often used, while the terms 
"active nematic" or "living liquid crystals" occur in the discussion of the orientationally ordered 
collective states of active particles [|39ll . 

Why study active systems? For one, such systems can display phases and phase transitions 
absent in systems in thermal equilibrium [|46l |47ll . For another, active matter often exhibits un- 
usual mechanical properties, including strong instabilities of ordered states to small fluctuations. 
Such instabilities may be generic in the sense that they should appear in any hydrodynamic theory 
which enforces momentum conservation and includes the lowest order contribution to the sys- 
tem stress tensor arising from activity Il44l . Fluctuations in active systems are generically large, 
often deviating qualitatively from the predictions of simple arguments based on the central limit 
theorem ll45l l46l l47l l37l . In common with a large number of related non-equilibrium systems, 
some examples of active matter appear to be self-tuned to the vicinity of a phase transition, where 
response can be anomalously large and fluctuations dominate the average behaviour. 




Figure 2: Cells of Escherichia coli labeled with a dye and examined in a fluorescence microscope. 
From these images, the behaviour of the flagellae of individual cells can be studied as E. coli 
swims ||49ll . Picture courtesy Prof. Howard Berg 

The study of active matter spans many scales. The smallest scales involve the modeling of 
individual motile organisms, such as fish or individual bacteria or even nanometer- scale motor 
proteins such as kinesin or dyneins which move in a directed manner along cytoskeletal filaments. 
Hydrodynamic descriptions of large numbers of motile organisms operate at a larger length scale, 
averaged over a number of such swimmers that is large enough for the granularity at the level of the 
individual swimmer to be neglected, while still allowing for the possibility of spatial fluctuations 
at an intermediate scale. 

Finally, all living matter is matter out of thermodynamic equilibrium. The source of this non- 
equilibrium is, typically, the hydrolysis of a NTP (nucleoside tri-phosphate, such as ATP or GTP) 
molecule into its di-phosphate form, releasing energy. This energy release can drive conforma- 
tional changes in a protein, leading to mechanical work. Biological matter is thus generically 
internally driven, precisely as demanded by our definition of active matter. Understanding the 
general principles which govern active matter systems could thus provide ideas, terminology and 
a consistent set of methods for the modeling of the mechanical behaviour of living, as opposed 
to dead, biological matter. This is perhaps the most substantial motivation for the study of active 
matter. 

The outline of this chapter is the following: Section[2]provides some examples of matter which 
can be classified as active, using our definition. Section [3] summarizes results on a simple model 
for an individual swimmer and summarizes necessary ingredients for a coarse-grained description 
of a large number of swimmers. Section |4]begins by illustrating the derivation of the equations of 



motion and the stress tensor of a simple fluid and then goes on to illustrate how this derivation may 
be extended to describe fluids with orientational order close to thermal equilibrium. It then goes 
on to discuss how novel contributions to the fluid stress tensor can arise from active fluctuations, 
deriving an equation of motion for small fluctuations above a pre-assumed nematic state. These 
calculations lead to predictions for the rheological behaviour of active systems and the demon- 
stration of the instability of the ordered nematic (or polar) state to small fluctuations. In Sect. [5j 
the theory of active gels is summarized briefly and some results from this formalism are outlined. 
Section [6] provides a brief summary. 

2 Active Matter Systems: Some Examples 

The subsections below list a (non-exhaustive) set of examples of active matter systems. The large- 
scale behaviour of some of these systems is the focus of later sections. 

2.1 Dynamical Behaviour in Bacterial Suspensions 

Dombrowski et.al. study the velocity fields induced by bacterial motion at the bottom of sessile and 
pendant drops containing the bacterium B. subtilis [fT2ll . These flows reflect the interplay between 
bacterial chemotaxis and buoyancy effects which together act to carry bioconvective plumes down 
a slanted meniscus, concentrating cells at the drop edge (sessile) or at the drop bottom (pendant). 
The motion exhibited by groups of bacteria as a result of this self-concentration is large-scale at 
the level of individual bacteria, exhibiting vortical structures and other complex patterns as a con- 
sequence of the hydrodynamic interaction between swimming bacteria. Thus, these experiments 
highlight the crucial role of hydrodynamics, coupled with the motion of individual active swim- 
mers, in generating self-organized, large-scale dynamical fluctuations in the surrounding fluid. 

2.2 Mixtures of Cytoskeletal Filaments and Molecular Motors 

The influential experiments of Nedelec and collaborators combine a limited number of cytoskeletal 
and motor elements - microtubules, kinesin complexes and ATP - in a two-dimensional geometry 
in vitro finding remarkable self-organized patterns evolving from initially disordered configura- 
tions Il32ll . These include individual asters and vortices as well as disordered arrangements of such 
structures, in addition to bundles and disordered states at varying values of motor densities. These 
are self-organized structures, large on the scale of the individual cytoskeletal filament and mo- 
tor, which require ATP (and thus are non-equilibrium) in order to form and be sustained. Several 
hydrodynamical approaches to this problem have been proposed, as in Ref. Il23l |4T1 0. These 
experiments and associated theory are summarized in Ref. [fT8l . 

2.3 Layers of Vibrated Granular Rods 

The experiments of Narayanan et.al. take elongated rod-like copper particles in a quasi-two- 
dimensional geometry, and vibrate them vertically in a shaker OTI . This agitated monolayer of 



particles is maintained out of equilibrium by the shaking, which effectively acts to convert vertical 
motion into horizontal motion via the tilting of the rod. The particles are back-front symmet- 
ric, and are thus nematic in character. These experiments see large, dynamical regions which 
appear to fluctuate coherently. Similar behaviour is predicted in theories of flocking behaviour 
induced purely by increasing the concentration in dense aggregates of particles held out of equi- 
librium |3T1 |48l . These experiments provide a particularly striking signature of non-equilibrium 
steady state behaviour in "flocking" systems: the presence of number fluctuations in such systems 
which scale anomalously with the size of the region being averaged, i.e., with a power of N which 
exceeds the central limit prediction of a \/N behaviour of number fluctuations [|48l . 

2.4 Fish Schools 

Experiments of Makris and collaborators use remote sensing methods on a continental scale to 
access the structure and dynamics of large-scale shoals of fish, containing an order of a few mil- 
lion individuals [25 J. Apart from the relatively rapid time-scale for the reorganization of the shoal 
- around 1-10 min in their experiments - these experiments also see evidence for "fish waves"; 
propagating internal disturbances within the shoal which occur at relatively regular intervals, rep- 
resenting disturbances at scales far larger than that of the individual fish. The speeds of such 
waves are larger, by around a factor of 10, than the velocities of the swimming fish. They appear 
to represent "locally interconnected compaction events" similar to the Mexican waves exhibited 
through the coordinated motion of spectators in stadia. It is interesting that similar wave-like exci- 
tations are predicted in the hydrodynamic theories of Refs. ||47ll44l . where they involve waves of 
concentration and splay [|47ll or splay-concentration and bend Il44l . 

2.5 Bird Flocks 

The STARFLAG collaboration has imaged large flocks of starlings (between 10 3 and 10 4 individ- 
uals at a time), using computer-aided imaging techniques to understand the dynamics of individual 
birds and how this dynamics is influenced by the spatial distribution and behaviour of neighbouring 
birds. Surprisingly, and contrary to what might have been naively expected, birds appear to adjust 
to the motion of the flock by measuring the behaviour of topologically (and not metrically) related 
neighbours O. Thus, at each instant, each bird appears to be comparing its instantaneous position 
and velocity to those of the 5-7 birds closest to it, making the adjustments required to maintain the 
coherence of the flock. This strategy appears to have the advantage that reducing the density of the 
flock should not then impact the coherence of the flock, since only topological and not metrical 
relationships are involved, a fact that calls into question the gradient expansions favoured by most 
theoretical work which represents flocking behaviour by coarse-grained equations of motion for a 
few hydrodynamic fields. 

2.6 Marching Behaviour of Ants and Locusts 

Ian Couzin's group at Princeton has investigated the transition to marching behaviour in locusts. 
Swarms of the desert locust, Schistocerca gregaria, in their non-flying form, can exist in a rela- 



tively solitary individualistic state as well as a gregarious collective state. In the "gregarious" state, 
such locusts can form huge collective marching units which forage all vegetation in their path. This 
has a huge social and economic impact on humans, affecting the livelihood of one in ten people on 
the planet in plague years. These experiments provide an experimental example of the collective 
transition in simple computational models of flocking behaviour. From more recent work from this 
group, it appears, somewhat unusually, that the transition is induced by "cannibalistic" behaviour, 
in which a locust is successful in biting the rear-quarters of the locust immediately in front 1171 [Toll. 

2.7 Listeria monocytogenes Motility 

The bacterium Listeria monocytogenes is a simple model system for cell motility, which derives 
its ability to move from the polymerization of actin, leading to the formation of a "comet tail" 
emerging from the rear of the bacterium Il3~5ll . Motility appears to arise from the deformation of the 
gel formed by the actin and cross-linking proteins as a consequence of continuing polymerization, 
which results in a propulsive force on the bacterium. The non-equilibrium comes from ATP-driven 
actin polymerization. Interestingly, many features of the experiment can be reproduced in in vitro 
systems where actin polymerization is initiated at the surface of specially treated beads, which then 
exhibit symmetry -breaking motility |f29l . 

2.8 Cell Crawling 

There is a vast and intriguing literature on cell crawling on substrates (81 HI Q31 [50J . Such crawl- 
ing appears to have four basic steps: the extension of cellular protrusions, the attachment to the 
substrate at the leading edge, the translocation of the cell body and the detachment at the rear. 
Cell crawling appears to be largely mediated by a meshwork of actin in gel form, whose fluidity 
is actively maintained through the action of myosin motors and other associated proteins ||8ll50ll. 
Interestingly, motility has also been observed in nucleus-lacking cell fragments excised by lasers, 
over a period of several hours. 

2.9 Active Membranes 

Experiments on fluctuating giant vesicles containing bacteriorhodopsin (BR) pumps reconstituted 
in a lipid bilayer indicate that the light-driven proton pumping activity of BR amplifies membrane 
shape fluctuations l|26l |27l . The BR pumps transfer protons in one direction across the membrane 
as they change conformation upon excitation by light of a specific wavelength. These experiments 
have been described in terms of a non-equilibrium "active" temperature ll26l 1271 . Hydrodynamic 
theories which interpret the experiments via a description of a membrane with a density of em- 
bedded diffusing dipolar force centres, i.e. as an "active membrane", have also been studied in 
some detail. It is interesting that several developments in the theoretical description of generic 
active matter reflect ideas first introduced in the context of active membranes |[33l [381 1421 . A de- 
tailed study of a biologically relevant model of an active membrane system which contains a large 
number of references to the literature is available in Ref. ETI . 



3 The Swimmer: Individual and Collective 



The examples provided in the previous section are indicative of some of the diversity exhibited by 
systems which fall under the general category of active matter. These can broadly be classified into 
systems in which the role of mechanical conservation laws (e.g. for momentum) are important and 
systems in which, typically, only the conservation of number is relevant. 

The flocking model of Vicsek and collaborators [j5TTl . recast in terms of hydrodynamic equa- 
tions of motion by Toner and Tu R6l l47l l48l . is an agent-based model in which cooperativity is 
driven by a density dependent interaction which tends to orient individual units in the direction 
in which their neighbours move. The environment - the ground and vegetation for the moving 
locusts, for example, provides merely the background in this case and has no other dynamical sig- 
nificance. In the case of the flock of starlings the communication between individual starlings does 
not appear to be primarily via the medium, but through visual contact. 

The case of the motile bacterium and the swimming fish, on the other hand, are cases where 
the medium plays an important role in transferring momentum to and from the swimmers and in 
determining the interactions between swimmers. This case will be discussed below, first for the 
situation of the individual swimmer and then for collections of swimmers, modeled via a hydrody- 
namic approach. 

3.1 The Individual Swimmer 

The central idea in the modeling of individual swimmers is that the system must be force-free, 
when averaged over lengthscales larger than the characteristic dimension of the swimmer 11391 1611. 
This is a consequence of Newton's third law, which imposes that the force exerted by the swimmer 
on the fluid must be equal and opposite to the force exerted by the fluid on the swimmer. Thus, if 
one considers the swimmer as a source of forces locally within the fluid which act upon the fluid, 
such a source cannot have a monopole component but may have a dipole (or higher multipole) 
component. 

Depending on the character of the swimmer, more distinctions are possible J6l[39l. Contractile 
swimmers or pullers (such as bacteria propelled by flagella at the head of the organism) pull fluid 
in along the long axis and push it out along an axis normal to their midpoint. Tensile swimmers 
or pushers push fluid out along their long axis and pull fluid in along the midpoints. They are 
propelled from the rear, justifying the terminology of pushers. 

Consider a simple model for a microscopic swimmer. In the Stokes limit, time does not enter 
these equations explicitly and the velocity field is completely specified by the boundary conditions 
imposed on the flow. (The use of the Stokes limit is generically justified in the case of bacteria, 
where characteristic Reynolds numbers at the scale of a single particle are of the order of 10~ 4 or 
smaller.) Reversing the velocity field at the boundaries should retrace the velocity field configu- 
ration, implying that the trajectory assumed by the swimmer in its configuration space cannot be 
time-reversal invariant. The helical or "corkscrew-like" motion of the flagellae of the bacterium E. 
coli (see Fig. [2]), discussed by Purcell 11341 . provides a particularly attractive example of how the 
limitations on directed motion at the low Reynolds numbers required for the Stokes limit approxi- 
mated to be valid can be overcome. 



Averaging over multiple strokes of the swimmer simplifies the description: the broken tempo- 
ral symmetry required for translation in a flow governed by the Stokes equations can be replaced 
by a broken spatial symmetry. This singles out the direction of motion of the swimmer. Following 
a model introduced recently by Baskaran and Marchetti (whose notations and treatment we follow 
closely in this section), the swimmer is modeled as an asymmetric, rigid dumbbell [6J. This dumb- 
bell consists of two differently sized spheres, of radii a L (large) and a s (small), forming the head 
and the tail of the swimmer. The length t of the swimmer is the length between the two centers. 
The orientation of the swimmer is given by the unit vector v, drawn from the smaller sphere to the 
larger sphere. 

Thus, the equations of motion of the two spheres, with locations r La and r# Q , are given by 

d t r La = u(r La ) 

d t v S a = u(r Sa ) , (1) 

where the constant length between the two spheres is imposed by the constraint that r^ Q — rsa = 
tv a and the velocity field is given by u(r). The no-slip condition at the surfaces of the spheres 
requires that the sphere move with the velocity of the fluid next to it. 

The velocity field obeys the Stokes equation and is constrained by incompressibility, 

?]V 2 u(r) = Vp + ^ active + f Noise, 

V-u(r) = 0, (2) 

where the "force" terms which enter on the right hand side are composed of a term which is associ- 
ated purely with activity (i.e. vanishes in thermal equilibrium) as well as of a second "noise" term 
modeling the fluctuations in fluid velocity arising from purely thermal fluctuations. Such noise 
terms do not conventionally enter the Navier-Stokes equations, but must generically be included in 
a coarse-grained description. Detailed expressions for these forces are available in Ref. J6]|. 

To solve the Stokes equations, we insert a delta function force 5(r)F on the right hand side, 
representing a fundamental source term from which more complex force configurations can be 
constructed. In Fourier space, the incompressibility condition is imposed as q • u(q) = 0, while 
the Stokes equations are — 7/g 2 u(q) + zqp(q) = F(q). This then gives 

u(r) ^f F _(aJ>^, (3) 

with 



q- F 

p(q) = i — (4) 



These are inverted by 



1/5 rr 

inr] \r r 

The 5 symbol is the unit tensor. In a Cartesian basis of i, j, k it is 



u(r) = — (- + -)• F. (5) 



8 = ii + jj + kk 



(6) 



The quantity acting on the force on the right hand side is the Oseen tensor, defined through 

O ij {r) = {5 ij +f i f j )/87c V r (7) 
for |r| > a^s, with f = r/|r| a unit vector. The Stokes equations are solved by the superposition 

Ui(r) = /y~] [Oij(r - r La ) - Oij(r - r Sa )] v a j ■ (8) 

a 

The divergence at short distances is eliminated through the definition: Cjj(|r| < ol.sO = 5ij/CL,s, 
where (l,s = Qnr]a LiS - 

The dynamics of an extended body in Stokes flow follows from translation of the hydrodynamic 
centre and rotations about the hydrodynamic centre. (The hydrodynamic centre refers to the point 
about which the net hydrodynamic torque on the body vanishes; it thus plays the same role in 
Stokes flow as the centre of mass in inertial dynamics.) For the problem of a rigid dumbbell in an 
external flow, the hydrodynamic centre is obtained from 

c + Csr s a L r L + a s r s 

r = 7 i 7- = T ' ( 9 ) 

(l + Cs a L + a s 

The equations of motion for the translation and rotation of the hydrodynamic centre follow from 

d t v c a = v v a + -=J2 F <xP + r <*( t ) 

to a = -^^r^ + rf(t), (10) 

where the angular velocity describing rotations about the hydrodynamic centre is defined by 

d t 9 a = 9 a x u a (11) 

and 

C = ((L + Cs)/2,(R = t 2 C- (12) 

The random forces r a and lead to diffusion at large length scales 10. 

The forces F a p and the torques r a p arise from hydrodynamic couplings between swimmers. 
An isolated swimmer is propelled at speed 

v ° = ~q — F (13) 

with velocity vq£>. This velocity arises purely as a consequence of the fact that the hydrodynamic 
and the geometric centers do not coincide. For symmetric swimmers, this velocity is zero and the 
swimmer is a "shaker" as opposed to a "mover". 

The interactions between swimmers can be calculated in the dilute limit by a multipole expan- 
sion, yielding 

F 12 « 2fa£ [3(r 12 • v 2 f - l] ^ (14) 

r l2 



for the hydrodynamic force exerted by the /3th swimmer on the ath one. In addition, expressions 
for the hydrodynamic torque between swimmers in the dilute limit can be derived and are presented 
inRef. HI. 

The hydrodynamic force decays as l/rf 2 , as follows from its dipole character. The torque 
consists of two terms: One is nonzero even for shakers and aligns swimmers regardless of their 
polarity. The second term vanishes for symmetric swimmers, serving to align swimmers of the 
same polarity. A more detailed discussion of the structure of the forces and torques for pushers 
and pullers is available in Ref. |[6]|. 



3.2 Multiple Swimmers 

The coarse-grained version of the many swimmer problem is defined through the following local 
fields First, we must have a density of active particles, defined microscopically in terms of 

c(r,t) = <£*(*-i£(*))>- (15) 

a 

Second, in the case in which the axes of a fore-aft asymmetric swimmer are largely aligned 
along a common direction - as in a magnet - we can define a local field describing polar order in 
the following way: 

P(r,*) = ^y<EM(r-r2(t))>. (16) 

Third, and finally, when one considers ensembles of interacting swimmers, we must also con- 
sider the possibility of additional and more subtle macroscopic variables representing orientational 
order. An ensemble of individual particles, each aligned, on average, along a common axis, is fa- 
miliar in soft condensed matter physics. Such systems are referred to as nematics and the ordering 
as nematic ordering. In such nematic order, the alignment is along a common axis but a vectorial 
direction is not picked out. 

Orientational order in the nematic phase is generally described by a second-rank, symmetric 
traceless tensor Q a p(x., t), defined in terms of the second moment of the microscopic orientational 
distribution function. This (order-parameter) tensor can be expanded as 

Qa(3 = 7^S {uaTlp - ^a/jj + (IJp ~ m a mp) . (17) 

The three principal axes of this tensor, obtained by diagonalizing Q a p in a local frame, specify 
the direction of nematic ordering n, the codirector 1 and the joint normal to these, labeled by 
m. The principal values S and T represent the strength of ordering in the direction of n and m, 
quantifying, respectively, the degree of uniaxial and biaxial nematic orders. 

In thermal equilibrium, the energetics of Q a p is calculated from a Ginzburg-Landau functional, 
first proposed by de Gennes, based on an expansion in rotationally invariant combinations of Q a p 



and its gradients IfTTI . The Ginzburg-Landau-de Gennes functional F is 



F = [ d^ATrQ 2 + \BTrQ 3 + jC(TrQ 2 ) 2 (18) 

J Z o 4 

+ E'iTrQ 3 ) 2 + h 1 (d a Q^)(d a Q^)]. 

Here, A = A (l — T/T*) T* denoting the supercooling transition temperature, A is a con- 
stant, L\ is an elastic constant and a, /3, 7 denote the Cartesian directions. Other elastic terms can 
also be included; this simple approximation corresponds to what is called the one Frank constant 
approximation. 

Two simplifications are possible and often convenient. First, we may assume uniaxial rather 
than biaxial order, since this is by far the more common form of ordering. In the ordered nematic 
state, the average orientation occurs along a direction n. This is the nematic director, defined to be 
a unit vector. We can thus work within a description in which the components of Q a/ 3 are written 
out in terms of the components of the nematic director n. This gives us 

3S ( 1 



Qaf3 = -y yiiUj - -Sijj . (19) 

The energetics of small deviations from the aligned state is obtained, in this representation, from a 
Frank free energy appropriate to uniaxial nematics: 

11 1 
f FO = -K^div n) 2 + -K 2 (n ■ curl n) 2 + -K 3 (n x curl n) 2 . (20) 

Here K\ is the splay elastic modulus, associated with a splay deformation V ■ n, K 2 is the twist 
elastic modulus and K 3 is the bend elastic modulus. The Frank constants K\, K 2 and K 3 have the 
dimensions of a force and can be represented as the ratio of an energy to a length scale. We can 
assume Ki ~ ksTc/a where a is a molecular length of order 1 nm. Note that this description uses 
three elastic constants (which can be reduced to the single elastic coefficient of Eq.[l9]by assuming 
that K 1 = K 2 = K 3 . 

With this background, the definition of a local field representing nematic order follows from [fTTTl . 



Qij = ^7^<I>-^ - ~ *£(*))> • (2D 



4 Hydrodynamic Approaches to Active Matter 

A class of questions relating to the modeling of active matter is concerned with (a) whether forms 
of collective ordering are at all possible in ensembles of interacting active particles and (b) whether 
such ordering, if assumed to preexist, can be shown to be stable against fluctuations. Our defini- 
tion of the single swimmer associated a direction with the swimming motion, the direction of the 
axis formed by connecting, say, the center of the small sphere to the center of the large sphere. 



The question is thus whether the axes or orientations of different swimmers can be aligned as a 
consequence of their interaction. 

In subsections below, the problem of deriving equations of motion for the nematic order param- 
eter field and the construction of a stress tensor appropriate to a nematic fluid are briefly examined. 
Results for the active nematic are summarized and the significance of these results for the Theolog- 
ical properties of active matter briefly outlined. 

4.1 Equations of Motion: Fluid and Nematic 

Identifying conservation laws and broken symmetries is the crucial first step in constructing hy- 
drodynamic equations of motion for the relevant fields in the problem. For both conserved and 
hydrodynamic fields, the relaxation of long wavelength fluctuations proceeds slowly, with the rele- 
vant timescales for the relaxation of the fluctuation diverging as the wavelength of the perturbation 
approaches infinity. 

For a simple fluid with no internal structure, the conservation laws for the local energy e, the 
density p and the three components of the momentum density gj are 



where j e is the energy current and it^ is the momentum current tensor, related to the stress tensor. 
The conserved momentum density itself acts as a current for another conserved density, the mass 
(equivalently, number) density, a relation which is responsible for sound waves in fluids. 

The hydrodynamic description of fluids with internal order (such as the nematic or polar fluid) 
must account for additional hydrodynamic modes arising out of the fact that the ordering represents 
a broken symmetry. For small deviations from equilibrium, one derives an equation for entropy 
generation and casts it in terms of the product of a flux and a force. Such fluxes must vanish at 
thermodynamic equilibrium. Close to equilibrium, it is reasonable to expect that fluxes should 
have a smooth expansion in terms of forces. 

As an illustrative example, consider the simple fluid in the absence of dissipation. We have, 
with u the velocity field, 



de 

dt 
dp 

dt 
dgj 
dt 



(22) 



g 



pu , 




pSij + ujgi = —Gij + pUiUj 



(23) 



The mass conservation equation is just 



dp 
dt 



V-(pu), 



(24) 



while the momentum conservation equation is 

dgi dpu, 



m m -VjTTij = -Vip-VjipuiUj) . (25) 
This is Euler's equation, usually written as 

^ + (u-V)u = — Vp. (26) 

at p 

The dissipative contribution to the stress tensor is accounted for by adding a term to the stress 
tensor, 

TTij = p5ij + pUiUj - a'^ . (27) 

Dissipation can only arise from velocity gradients, since any constant term added to the velocity 
can be removed via a Galilean transformation. The dissipative coefficient coupling the stress tensor 
to the velocity gradient is most generally a fourth rank tensor, 

°'%3 = Vijkl^kUl . (28) 

However, symmetry requires that 

2 

a ij = Vv^iUj + VjUi - -SijV ■ u) + C<%V • u . (29) 
This gives the Navier-Stokes equations. Assuming incompressibility, we have 



d ± = = -V • (pu) = -V • u . (30) 

Thus, 

p^ + p(u- V)u = -Vp + r?V 2 u, (31) 

along with the constraint V ■ u = 0. The velocity field thus has purely transverse components. 

How is this to be generalized for a nematic fluid? First, for a nematic fluid, we must have 
an equation of motion for the director n (or, equivalently, for the Q Q( g tensor) in addition to the 
equations for the conservation of matter, momentum and energy IfTTl |43l . Second, distortions 
of configurations of the nematic order parameter field also contribute to the stress tensor of the 
system HHIH. 

The director is aligned with the local molecular field in equilibrium; local distortions away from 
the molecular field direction must relax in order to minimize the free energy. The local molecular 
field is defined in the equal Frank constant approximation as 

hi = KV 2 rn . (32) 

Also, the director does not change under rigid translations at constant velocity. Thus, the leading 
coupling of n to u must involve gradients of n. This can then be written as 



dt 



-^r - XijkVjU k + X i = , (33) 



where X' is the dissipative part of the current. This dissipative part can be written as 



K = , (34) 



where 7 is a dissipative coefficient and the projector isolates components of the fluctuation in the 
plane perpendicular to the molecular field direction. 

The constraint n ■ dn/dt = implies that there are only two independent components of the 
tensor Xij k - These can be taken to be symmetric and antisymmetric, defining 

Xijk = ^A(5jn fc + tfarij) + -A 2 (5jn fe - S^rij) , (35) 

where 

Sfj = dij - riiUj . (36) 

Under a rigid rotation, 

On ■ 1 

— - = w x n = -(V x u) xn, (37) 

mandating that the coefficient A 2 of the antisymmetric part must be -1. 

Thus the final equation of motion for the director, the Oseen equation, takes the form [fTTll43l 



d t rii + v ■ Vrii + UijUj = Sf, ( Xu^nj + — ) . (38) 



hi 
7 

where Uij = (1/2) (diUj + djUi). 

The stress tensor of the nematic consists of three parts. The first is the thermodynamic pressure 
p, while the second is the viscous stress, given by the tensor 

cr' i:j = aininjn k niA k i + a 2 n j N i + a 3 niNj + a^A tj + a 5 njn k A ik + a 6 nin k A jk , (39) 

constructed from symmetry allowed terms, where A represents the symmetric part of the velocity 
gradient as before and the vector iVj = hi + |(n x curluji is the change of director with respect 
to the background fluid. One relation, due to Parodi, connects the coefficients a\ . . . a§; there are 
thus 5 independent coefficients of viscosity in the nematic ffTTTl . 

The third contribution to the stress tensor is the static (elastic) contribution arising from defor- 
mations in the director field [1 1 J, i.e., 

ff « = -^ v <"< ■ (40) 



A more detailed discussion of the equation of motion of the nematic order parameter and the stress 
tensor is available in Ref. ifTTI . 



4.2 Active Orientational Order and its Instabilities 



The central idea behind the modeling of the active nematic is that out of thermal equilibrium, new 
terms enter the equation of motion for the nematic director as well as the stress tensorll] These 
terms are more "relevant" than the terms mandated by thermodynamic approaches, in the sense 
that their effects are stronger at long wavelengths and at long times, i.e. in the thermodynamic 
limit. To demonstrate this, assume nematic or polar ordered particles, given by a unit director field 
n. The "slow" or hydrodynamic variables are (i) the concentration fluctuations 5c(r,t), (ii) total 
(solute and solvent) momentum density g(r, t) = pu(r, t) and broken symmetry variables whose 
fluctuations 5n± — n — z. For polar ordered suspensions, there is a non-zero drift velocity v n. 
The momentum density evolves via 

dgj 
dt 

with (Jij the stress tensor. 

What is special about self-propelled particle systems is that the active contribution to the stress 
tensor is proportional to the nematic order parameter, i.e., 



-VjtTii , (41) 



afj oc ( riiUj - ^5ij ) . (42) 



Pascals law is thus violated in the active nematic, since it has a non-zero deviatoric stress. Also, 
crucially, terms obtained from the standard near-equilibrium analysis are of higher order, since 
they involve gradients of the nematic order parameter. 

This result can be derived from a microscopic calculation using the simple model for a single 
swimmer discussed above, by simply Taylor expanding the dipole term which appears in the force 
density, about the hydrodynamic centre. This yields 

°$j = ^4^fc(r, t) L inj - , (43) 

where ol and a# were defined earlier in the context of the single swimmer. We assume an initially 
aligned state in which the director points, on average, along the z direction. Fluctuations away 
from this are given by Sn± = n — z. Fluctuations about the averaged value can be parameterized 
by 




5n = y 5n y \ . (44) 

Also, we can expand the velocity field u in terms of fluctuations about its mean value u , 

Su x 

riu = | Su y J . (45) 

u + 5u z 



'in this section, we follow the treatments of Refs. i44l[T51 closely. 



Finally, we must allow for concentration fluctuations, via 

c(r, t) = c + Sc(r, t) . (46) 
Using the equation of motion for the velocity field, defined through u = g/p, we get 



(a L + 05) / 1 

;-)(c + ()c(r,i))(n i n j - -%) 



2 V v '" v * J 3 



(47) 



This then gives, once we insert the expressions for small fluctuations and retain terms to the lowest 
order 

dtSu x = (yClL + — {-)c d z 5n x 
z p 



(a L + a s ) f 

o t 0Uy = {-)c d z bn y 

Z p 



d t 5u z = K ± as) (t)c (djn x + d y 5n y ) + ( ° L ± as) { f -)d z 5c{v, t) , (48) 

Zi p Zi p 

which we can write as 

d t 5u± ~ w d z 5n± 

d t Su z ~ w V± ■ 5n± + ad z 5c(r, t) , (49) 

where a ~ fa/p and w ~ c a and a ~ ( a ^+ a s) ^ ^ e mus t now impose incompressibility, which 
requires that 

V±5u ± + d z u z = . (50) 

In Fourier space, this is 

iq± ■ Su ± (q,t) +iq z u z (q z ,t) = . (51) 

Incompressibility implies that the velocity field must be purely transverse, a condition that is easily 
imposed by using the transverse projection operator 

Pl f = (** - q ~f) ■ (52) 

Operating this projection operator on a vector field isolates those parts of the field which have no 
longitudinal component. Doing this yields 

2 

d t Su x = -iw (l - ^)q z 5n x - ^-^-q z bn y - ia^^(q x 5n x + q y 5n y + q z 5c) 

gZ gZ gZ 

2 

d t Su y = -iw (- ( ^-)q z 5n x + (1 \)q z bn y 

gZ gZ 

- ia^^(q x 5n x + q y 5n y + q z 5c) , (53) 



which can be written in the compact form 



= -iw q z (l - 2q_Lq_L/g 2 ) 5n ± - i-^a(q ± 5c) . (54) 
Number conservation implies 

85c 

-Vj, (55) 



dt 



where j = cv n = cv (5n± + z), so 



d d\ 

dt + V °dz) SC + C ° v °^- L ' Sn = • ^ 

The equation of motion for polar-ordered particles whose orientation is described by n contains 
a term representing advection by a mean drift, a term describing the consequences of a non- 
equilibrium osmotic pressure and other terms familiar from our brief description of neurody- 
namics close to equilibrium in the previous section lPRl . 

d5n ± d5n± _ . l 0/ <5u^ _ . 1 .d5u± _ . 

= - AlK »— - + 2 d{ i>r ~ + 2^ ( — + Vi " J • <57) 

Taking the curl of this equation as well as the equation for for gives coupled equations 
for the dynamics of twist V x and/or vorticity Vxui. Related results follow from taking the 
divergence of these two equations, resulting in coupled equations for V ■ ni, V • ui and V5c, which 
are complicated but have also been analysed p4|. These coupled equations have been shown to 
possess wave-like solutions. 

The results of these calculations are summarized in the following: 

• A linearized treatment, ignoring viscosity, for the polar or apolar cases at lowest order in 
wave number, yields propagating modes with a characteristic instability in the case of purely 
apolar active particles. 

• Retaining viscosity, in the steady Stokesian limit where accelerations are ignored, polar and 
nematic orders at small wave numbers are generically destabilized by a coupling of splay (for 
contractile particles) or bend (for tensile particles) modes to the hydrodynamic velocity field. 
In this limit, for nematic order, this instability has been referred to as a generic instability. 
(A physical picture for the origins of this instability is provided in Ref. 11391 .) 

• Any mechanism for introducing screening into the hydrodynamic interaction will suppress 
the instability, such as introducing boundaries into the system (6l[39l| or applying a shear 11301 . 



Number fluctuations in ordered collections of self-propelled particles are anomalously large [|46 
|47ll . The variance ((5N) 2 ), scaled by the mean N, diverges as iV 2 / 3 in three dimensions in 
the linearized treatment of Ref. 11441 . Physically, this is a consequence of the fact that (a) 
orientational fluctuations (director distortions) produce mass flow and (b) such fluctuations 
are large because director fluctuations arise from a broken symmetry mode. 



In a more general context, activity provides new currents both for matter and momentum, 
beyond those which would be predicted from a theory based on small perturbations away from 
thermodynamic equilibrium. It is these new currents which are the source of the novel results 
indicated above. 



4.3 Rheological Predictions 

The following predictions for the rheological properties of active nematics are obtained from |fT51 
and the discussion here follows that reference closely; see also lf24l [T[|22l. The active contribution 
to force densities within the fluid follows from 

F a = V • a a . (58) 

The activity contributes to traceless symmetric (deviatoric) stress through 

a a - \Tra a ~ WiQ + W 2 Q 2 , (59) 

where W% and W 2 reflect the strength of the elementary dipoles. The full stress tensor has contri- 
butions from the fluid (cr v ) as well as the order parameter field (cr OP ), giving 

a- = (T a + a v + a OP . (60) 

The time-dependence of a is assumed to be slaved to the time-dependence of the order parameter 
field Q. The equation of motion for Q, upon linearization, should take the form 

?Q. = _Iq + DV 2 Q + A A + . . . , (61) 

at t 

where r is an activity correlation time, D is a diffusivity, A is a kinetic coefficient and higher 
order terms have been dropped in favour of the lowest-order ones. Using these, we can calculate 
the stress response to a shear in the xy plane. In Fourier space, this is 



, (a + W)X 

—lUJ + T 1 

G\u) - iG"(u) 



A 



iA xv , (62) 



defining the storage and loss moduli G (u) and G (u). The important results which follow from 
this analysis are the active enhancement or reduction r] act oc Wr (depending on the sign of the 
parameter W) of the effective viscosity at zero shear rate. The theory also predicts strong vis- 
coelasticity as r increases. For passive systems W = 0. For active systems, W is non-zero and the 
storage modulus then behaves as 

G'(ut > 1) ~ W . (63) 

All these effects are expected to be enhanced if the transition is to a polar ordered phase, rather 
than to a nematic. 



5 Active Gels: Summary 



A parallel line of activity, centered on models for specific biological phenomena such as the 
symmetry-breaking motility exhibited by beads coated with polymerizing actin and the dynamical 
topological defect structures obtained in mixtures of motors and microtubules, outlines a descrip- 
tion of active matter in terms of active gels lfT3l[T7l[T9ll20l[T6ll . 

The philosophy of these approaches is the following: Rather than begin from a microscopic 
model for a swimmer or individual moving particle and then generalize from the microscopies to 
realize symmetry-allowed equations of motion for the fluid velocity field and for the local con- 
centration of swimmers, one can start with a coarse-grained continuum model for a physical vis- 
coelastic gel which is driven by internally generated, non-equilibrium sources of energy IfToll . The 
equations of motion for the stresses in this gel as well as for local order-parameter-like quantities 
describing, for example, polar order at a coarse-grained scale are constructed using the basic sym- 
metries of the problem. The non-equilibrium character of this problem follows from the fact that 
such equations of motion do not derive from an underlying free energy^. 

The passive gel has a viscoelasticity whose simplest representation is via a Maxwell model, 
exhibiting solid-like behaviour at short times and fluid-like behaviour at long times. In this model, 
the deviatoric stress a a p is related to the strain rate tensor u a p = \ [d a up + dpu a ), where u is the 
velocity field in the gel, via 

d -^ + ^=2Eu aP , (64) 
at r 

where E is a shear modulus obtained at short times. The simple time derivative must be augmented 
by convective terms, as well as terms representing the effects of local rotation of the fluid, to 
enforce Galilean invariance and the appropriate frame independence. 

Polar order in such gels is assumed to be weighted by a free-energy-like expression obtained 
from the theory of polar nematic liquid crystals. This takes the form: 



7 = \ dr 



-f (V ■ P) 2 + -y (P • (V x p)) 2 + -f(p x (V x p)) 2 + W-p-jp 2 



(65) 



There are three Frank constants for splay, twist and bend, as in the nematic case. The (non- 
zero) constant k is permitted by the vectorial symmetry of the polar case. The amplitude of local 
order is parametrized by the constant h?<. 

The hydrodynamic theory of active gels begins by identifying, along classical lines, fluxes 
and forces. The hydrodynamic description contains phenomenological parameters, called Onsager 
coefficients. These fluxes are the mechanical stress a a p associated with the mechanical behaviour 
of the cell, the rate of change of polar order (the polarization) P and the rate of consumption 
of ATP per unit volume r. The generalized force conjugate to the ATP consumption rate is the 
chemical potential difference Afi between ATP and the products of ATP hydrolysis, while the 
force conjugate to changes in the polarization is the local field h, obtained from the functional 
derivative of the free energy, i.e. h = —5J r /5p. The force conjugate to the stress tensor is, as 



2 A summary of these results can be found in [17], from where most of this material is drawn. 



usual, the velocity gradient tensor d a up. This can, as is conventionally done, be expanded into its 
traceless symmetric, pure trace and antisymmetric parts. A similar expansion can be made for the 
stress tensor. 

The next step is to construct equations of motion for the deviatoric stress, using the convected 
Maxwell model with a single viscoelastic relaxation time. The equation must couple the mechani- 
cal stress and the polarization field as well as include a term coupling activity to the stress. It takes 
the form 

2rjU a p = (l + T j^J jc-a/3 + CA/i<? a/ 3 + TA aj3 - ^(pah^ + ppha - -/l-yP-y^) j , (66) 

where the co-rotational derivative is 

D ( 9 9 \ 

jx**? = + M 7^r J a ^ + Kt^ + w /37 a 7°J ; ( 6? ) 

the tensor A a p describes geometrical nonlinearities arising out of generalizations of the Maxwell 
model to viscoelastic fluids, and q a p = \{j> a p$ — kp 2 8 a p)- The antisymmetric part of the stress 
tensor leads to torques on the fluid and is obtainable from 

The viscoelastic relaxation time is r, the coefficient v\ describes the coupling between mechanical 
stresses and the polarization field, while the parameter ( is the coefficient of active stress gener- 
ation, acting to couple activity to the stress. The second flux, defined from the rate of change of 
polarization is given by 

DP 

The Onsager relation for the polarization is 

^-p a = —h a + \ip a Afl - VtUapPp - PiUppPa . (70) 

Dt 7i 

These include several phenomenological parameters, such as the rotational viscosity 7 X which 
characterizes dissipation from the rotation of the polarization as well as the constants V\ and V\. 
Then, we must have an equation for the rate of consumption of ATP. This takes the form 

r = AA/i. + CPaP/3 u a/3 + C u aa + ^lPah a ■ (71) 

These are simple but generic equations representing the basic symmetries of the problem lfToTl . 
They can be shown to have interesting and surprising consequences: an active polar gel can exhibit 
spontaneous motion as a consequence of a gradient in the polar order parameter, as well as defects 
in the polar ordering which are dynamic in character [fTTl . These ideas have been applied to the 
study of the motion of the cell lamellipodium and to the organization of microtubules by molecular 
motors lfT6ll . 



The generality of these equations follows from the fact that they are motivated principally 
by symmetry considerations. Thus, even though they describe intrinsically non-equilibrium and 
highly nonlinear phenomena, for which the rules of constructing effective, coarse-grained equa- 
tions of motion for the basic fields are not as well developed as the theory for the relaxation of small 
perturbations about thermal equilibrium, the "unreasonable effectiveness of hydrodynamics" may 
well hold in their favour. 

6 Conclusions 

This chapter has provided a brief review of the field of what is currently called active matter. The 
emphasis has been on attempting to clarify the basic ideas which have motivated the development 
of this field, rather than the details of the often intricate and complex calculations implement- 
ing these ideas. Much recent and important work, including numerical calculations - illustrative 
references are Refs. [|40ll28l[36l - has been omitted entirely for the sake of compactness. 

As indicated in the introduction, the importance of this field would appear to be that it might 
suggest ways of thinking about the response and dynamics of living systems, while providing a 
largely self-consistent framework for calculations. Several non-trivial insights have already been 
obtained from these calculations, particularly in the identification of the generic instability of po- 
lar or nematically ordered states in the presence of the long-ranged hydrodynamic interaction, 
the connection between microscopic models and their hydrodynamic limits, as well as a compre- 
hensive theory of active gels, generalizing ideas from nematic physics. The precise relationship 
between macroscopic, symmetry-based hydrodynamic equations representing active nematics and 
an underlying microscopic theory has been substantially clarified, as in the work of Ref. [6] and 
references cited therein. To what extent further developments in this field may aid the increasingly 
active dialogue between physics and the engineering sciences on the one hand and the biological 
sciences on the other remains to be seen. 
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